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Abstract 

We investigate the fate of a relativistic star cluster with a dense core which 
is undergoing a gravothermal catastrophe and is far from thermodynamic 
equilibrium. Nonlinear cooperative contributions are included in the standard 
transport equations for the last stage of evolution of a highly dense core of 
stellar remnants. We find that the core redshift does not necessarily increase 
without limit as the core becomes increasingly dense, preventing collapse to 
a black hole. In particular, the redshift can remain less than the critical 
value for relativistic collapse, resulting in a stable, massive dark core with a 
Newtonian mantle and halo. 

Subject headings: black hole physics - galaxies: nuclei - gravitation 
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I. INTRODUCTION 



According to numerical integrations of the full Einstein equations for the dynamical evo- 
lution of an arbitrary spherical, collisionless system in general relativity (Shapiro and Teukol- 
sky 1985a,b,c, 1986), star clusters become relativistically unstable at sufficiently high central 
redshift, z c > 0.5, confirming the original speculation by Zel'dovitch and Podurets (1965) 
and perturbative calculations performed by Ipser and Thorne (1968), Ipser (1969,1970,1980), 
Bisnovatyi-Kogan and Thorne (1970), Fackerell (1970) and Fackerell et al., (1969). How- 
ever, more recent work by Rasio, Shapiro and Teukolsky (1988,1989), Bisnovatyi-Kogan, et 
al., (1993) and Merafina and Ruffini (1995), demonstrated that a family of cluster evolution 
models were relativistically stable for z c — > oo. Merafina and Ruffini proposed three possible 
configurations of cluster equilibrium models in terms of central redshift versus temperature 
diagrams. They found that the stable relativistic equilibrium models, corresponding to the 
possibility of reaching infinite central redshift, were not of astrophysical interest, because 
they have ratios M core /M ~ 10~ 24 and R core /R ~ 10~ 25 , where M core and -R cor e denote 
the core values of the cluster mass and radius, respectively. We shall accept the point of 
view that astrophysically interesting clusters, which can evolve dense relativistic cores with 
masses 10 7 — 10 9 M , will become unstable for z c > 0.5 and collapse to black holes. 

The question remains: do supermassive black holes form at the center of galaxies? In 
a recent article (Moffat 1997), a model of a galaxy was constructed from a relativistic 
massive core and a Newtonian mantle and halo. The core consisted of stellar remnants 
such as neutron stars, brown dwarfs or planets with individual masses, m ~ 1M . This 
would not contradict the observations obtained by Lugger et al., (1992), which ruled out 
a change in mass-to-light ratio from ~ 2 at 1 arsec to > 20 at 0.1 arcsec, yielding large 
broadband color gradients. Using the high-velocity maser emission data obtained for NGC 
4258 (M ~ 4 x 10 7 M Q )(Miyoshi et al, 1995; Maoz 1995), which impose constraints on the 
mass distribution of the dense stellar cluster, the evaporation time obtained was t cvap ~ 10 
Gyr, a time scale big enough compared to the Hubble time. 
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Calculations have been performed of dense (> 10 8 stars pc~ 3 ) galactic nuclei composed 
initially of main-sequence stars (Gold, Axford, and Ray 1965; Spitzer and Stone 1967; Col- 
gate 1967; Sanders 1970; Spitzer 1971). It was found that the evolution is mainly determined 
by coalescence and stellar collisions, and massive stars would undergo supernova explosions 
that can leave behind compact stellar remnants. The endpoint of this evolutionary phase 
is likely to be a cluster of neutron stars or other stellar remnants. Such a cluster could 
initially consist of 10 s compact objects, within a region < 0.1 pc in radius and moving with 
a velocity dispersion 800 — 2000 km s" 1 . Two epochs characterize the subsequent Newtonian 
evolution of the cluster of compact stars, (i) a low-redshift (z c < z co \\ ~ 10~ 2 ) point-mass 
epoch during which the core undergoes collapse driven by the "gravothermal catastrophe" 
(Lynden-Bell and Wood 1968; Cohn 1980), (ii) a subsequent short, high-redhift (z c > z co \\) 
epoch dominated by coalescence and collisions of compact stars in the cluster. 



II. BASIC NONLINEAR TRANSPORT EQUATIONS 

We must now investigate the evolution of the dense core following these two epochs. 
Shapiro and Teukolsky (1987), Goodman and Lee (1989), Quinlan and Shapiro (1989), Lee 
(1995), Quinlan (1996), van der Marel et al., (1997) have studied this epoch using transport 
equations for a cluster of compact stars (Lightman and Shapiro 1978). The number of core 
stars decreases and, at the same time, the central velocity dispersion and redshift increase 
ultimately leading to a relativistic redshift, z c > z cv \ t ~ 0.5, and therefore to gravitational 
collapse to a black hole. 

These authors assumed that the transport equations were linear differential equations 
in the time variable t. The arguments against dark clusters of compact stars rely mainly 
on the very high density of the core (> 10 8 M o pc~ 3 ). Treated as a fluid, the gravothermal 
catastrophe epoch is expected to be far from thermodynamic equilibrium (Lynden-Bell and 
Wood 1968). 

There are many examples in physics in which subsystems cooperate with each other in a 
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self-organizing way (Haken 1975; Haken 1983; Bak, Tang and Wiesenfeld 1988; Bak 1996). 
The behavior of the total system can show characteristic changes which can be described 
as a transition from disorder to order, or a transition from one state to another. There are 
numerous examples of such behavior at both the macroscopic and quantum levels. Pro- 
nounced cooperative phenomena may occur in physical systems far from thermodynamic 
equilibrium. Such physical systems far from equilibrium can display ordered states created 
and maintained by an energy flux passing through the system. The choices of order param- 
eters, well-known in phase-transition theory, play an important role. The order parameters 
represent the behavior of the system on a macroscopic scale and therefore describe macro- 
scopic variables. Normally the order parameters satisfy simple differential equations with 
respect to the time variable, for the relaxation time of order parameters is normally much 
greater than those of the subsystems. The variables of the subsystem can be eliminated 
without increasing the degree of the time derivatives. 

We shall generalize the differential equations for the dense core given by Lightman and 
Shapiro (1978) and Shapiro and Teukolsky (1985) to the following equations: 

E c = -( ai /t r )E c - (Pi/t 3 r )E 3 c + (l/t coll )E c , (la) 
N c = -(a 2 /t r )N c - {f3 2 /tl)Nl - (l/t coll )N c , (lb) 
7h c = (l/t co ii)m c , (lc) 

where N c (t) is the number of core stars, N c = dN c /dt, E c (t) is the characteristic core binding 
energy and m c (t) is the mean stellar mass. By using the virial theorem we get 

E c « hs c m c vl « ^{N c m c ) 2 /r c , (2) 

and the relations 

n c = iV c /(47rr c 3 /3) « 2 x 10 V -3 ^^^. 3 ). (3a) 

r c = 0.2 pc (N c ^v~lm c ^), (3b) 

where N c>8 = N c /10 8 ,v Ct s = f c /10 3 km s _1 , and m c ^ = m c /M & . The central relaxation time 
is (Spitzer and Hart (1971)): 
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3\ 1/2 

4nm 2 c n c ln(0.4iV c ) 



0.8xl0Vr(Ar>-|m c ,A" 1 ) ) (4) 



,2, 

where A§ = ln(0.4AT c )/ ln(0.4 x 10 8 ). Moreover, the dynamical time scale td -C t r is 



t d « r c /w c « 200 yr (N cfi vJm Ci *), (5) 

and the two-body collision time scale is 

tcoii « « 6 x 10 13 yr (A^V^V,*), (6) 

where v oo ~ 2 1 / 2 t> c is the asymptotic relative velocity, and <7 co ii is the collision cross section. 

The physical source of the nonlinear contributions in the transport equations for the 
order parameters E c and N c could be energy pumped into the core from binary gravitational 
radiation, heat dissipative effects caused e.g., by tidal interactions that can transfer large 
amounts of energy to the core in its final stage of evolution, or post-Newtonian contributions 
in the self-gravitating interactions of the point masses in the core. 



III. APPROXIMATE SOLUTION FOR THE REDSHIFT 

By using the result for the redshift (Cohn 1980; Stuart and Teukolsky 1985): 

z c w C w Uv 2 c /3 «5x 10~ 5 ^ 3 , (7) 
where C is the central potential, we can replace Eq. ([Ta|) by 

z c = ( ai /t r )z c - (h/t 3 r )z 3 c + ((3 2 /t 3 r )N 2 z c + (l/t co ii)^, (8) 

where a\ = a 2 — a\, b\ = (3N c m c /28) 2 f3 1 and the total core mass N c m c is constant. 

We solve Eqs.(|), ( [lb]) and ( [Tc[ ) subject to the constraint: r < H^ 1 = 2 x 10 10 yr, where 
t is the remaining time for the gravothermal catastrophe to drive secular core collapse to 
completion (A^ — >• 0). A solution for aii and a 2 (Cohn 1980; Stuart and Teukolsky 1985), 
using a matching to Fokker-Planck calculations for the advanced core collapse in the absence 
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of collisions, yielded the result a\jai = 0.701 (ai > 0). It should be noted that the Fokker- 
Planck approximations used by Cohn (1980) eventually break down as the number of stellar 
remnants in the core decreases and large angle scattering and binary formation become 
important, so these caculations may not include nonlinear cooperative phenomena in the 
last highly dense stage of core evolution. 

We shall assume that the parameter ai remains positive for the nonlinear equations, 
Eqs.((|) and fllbl) , that t co n = oo, that the parameter (3 2 ~ and that in the time scale of 
the gravothermal catastrophe evolution t r ~ const. Then, Eq.(|) can be written as 



Z r . 



az c - bz 3 c , (9) 



where a = a\jt r and b = bi/t 3 r . This equation has an exact time dependent solution for 
a > and b > given by 

a 1 / 2 



c {b + Caexp[-2a(t-t')}} 1/2 ' 
where C is an integration constant. For b = we get 



(10) 



z c = C- 1/2 exp[a(t-t')], (11) 

and z c — > oo as t — > oo, which leads to the catastrophic collapse to a black hole. However, 
for b > we get as t — ► oo: 

z c w (a/6) 1 / 2 , (12) 

and if (a/b) 1 ^ 2 < 0.5, then the dark cluster can be relativistically stable if the adiabatic 
index I\ w 5/3 and the cluster has a Newtonian mantle and halo with a mantle radius, 
r m > 0.1 pc (Moffat 1997); the massive core will contain most of the mass of the galaxy. 



IV. DISCUSSION 



Observations of nuclear regions in the centers of galaxies have detected high velocities of 
revolution, and luminosity spikes suggesting masses between 10 7 and several 10 9 M & (Ford 
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et al., 1994; Miyoshi et al., 1995; Kormendy and Richstone, 1995; Kormendy et al., 1997). 
However, masses for more than 60 other nearby galaxies are of order ~ 10 6 M & (Kormendy 
1987). For our Galactic center, the mass may be significantly less than 1O 6 M (Allen and 
Sanders 1986; Kundt 1990; Sanders 1992). The black hole models of clusters have been 
used to explain AGNs, and other intense extragalactic radio sources powered by jets of 
plasma and magnetic fields produced by the compact massive black hole at the center of 
the galaxy (Begelman, Blanford, and Rees 1984). However, Kundt (1996) has formulated 
an alternative model to explain AGN observations, based on an application of standard 
accretion disk theory to the central parts of the galactic disk. A burning disk can possibly 
explain the phenomena commonly attributed to a supermassive black hole. It is not clear 
at present whether these models can successfully describe the complex phenomena of AGNs 
for all the observed galactic nuclei. 

We have demonstrated that by using nonlinear transport equations for the order param- 
eters E c and N c , which take into account nonlinear cooperative effects of the subsystem 
of point masses, the inner core of a cluster may become relativistically stable, after hav- 
ing been formed by an epoch of gravothermal catastrophy. Treated as a fluid, the dense 
core is far from thermodynamic equilibrium and from known physical processes, we can 
expect cooperative phenomena to be important. The nonlinear cooperative contributions 
to the differential equation for the redshift can prevent it from achieving the critical value, 
z c ~ ^crit ~ 0.5. In the linearized approximation to the transport equations, the standard 
result follows that as gravothermal catastrophe develops, the density of the core increases 
and the number of core stellar remnants decreases until z c > 0.5, at which point the core 
inevitably undergoes gravitational collapse to a black hole. 

We have used a heuristic model of the nonlinear cooperative effects in a highly dense core 
far from thermodynamic equilibrium, in order to obtain a picture of the qualitative behavior 
of the last epoch of core evolution. A more complete solution would require a numerical 
analysis of the problem. Equations such as Eq.flTaD should include stochastic contributions, 
leading to non-zero correlation functions between two-particle interactions, and use should be 
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made of a Fokker-Planck equation and a related probability density distribution to analyze 
the dynamical evolution of the relativistic core. Moreover, it is unrealistic to ignore collisions, 
binary formation and gravitational radiation (Quinlan and Shapiro 1989) in the last stage 
of core evolution. 

In spite of the shortcomings of our qualitative picture of the last stage of core evolution, it 
is clear that non-linear self-organizing, cooperative effects, associated with phase transition 
phenomena for the highly dense cluster cores, can dramatically change the collapse scenarios 
and possibly prevent the formation of supermassive black holes. 
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